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CANONICAL FORMS FOR NONLINEAR SYSTEMS 
Renjeng Su,* L. R. Hunt, + and George Meyer 
Ames Research Center 
SUMMARY 

Necessary and sufficient conditions for transforming a nonlinear system to a 
controllable linear system have been established, and this theory has been applied 
bo the automatic flight control of aircraft* These transformation results show that 
the nonlinearities in a system are often not intrinsic, but are the result of unfor- 
tunate choices of coordinates in both state and control variables. Given a nonlinear 
system (that may not be transformable to a linear system) , we construct a canonical 
form in which much of the nonlinearity is removed from the system. If a system is 
not transformable to a linear one, then the obstructions to the transformation are 
obvious in the canonical form. If the system can be transformed (it is called a 
linear equivalent), then the canonical form is a usual one for a controllable linear 
system. Thus our theory of canonical forms generalizes the earlier transformation 
(to linear systems) results. Our canonical form is not unique, except up to solu- 
tions of certain partial differential equations we discuss. In fact, the important 
aspect of this paper is the constructive procedure we introduce to reach the canoni- 
cal form. As is the case in many areas of mathematics, it is often easier to work 
with the canonical forms than in arbitrary coordinate variables. 


I. INTRODUCTION 

Suppose we have a nonlinear system 

x(t) = f (x(t) ) + u (t) g (x(t) ) (1) 

where f and g are real-analytic vector fields on 3R n and f (0) = 0. If 
f(x) = Ax and g(x) = b, where A is nxn and b is nxl, and b,Ab, . . . ,A n-1 b 
are linearly independent, then we can always find new coordinates y 1 ,y 2 , • • * »Yn» v > 
where the y ± are functions of x, and v is a function of (x,u), such that the 
system becomes a "string of integrators" 

Yi - y 2 

y 2 - y 3 

( 2 ) 

^n-1 = y n 

y n - v 


*Research Associate of National Research Council at NASA Ames Research Center. 
^Research supported by NASA Ames Research Center under the IPA Program and the 
Joint Services Electronics Program at Texas Tech University, Lubbock, Texas, under ONR 
Contract N00014-76-C-1136. 



This is a canonical form for the controllable linear system. 

We let [f ,g] , (ad 2 f ,g) , . . ., (ad k f,g), . . . denote successive Lie brackets 
(definitions are given in the second section) of the vector fields f and g from 
equation (1). If g,.[f,g], . . . , (ad n-1 f , g) are linearly independent and 

g> [f jg ] (ad n-2 f,g) are involutive, then from reference 1 we know there exists 

a neighborhood of the origin in x-space and new coordinates y x ,y 2 > • • • >y n and 
control v so that we have equation (2) . Thus the system (1) appears as a nonlinear 
system only because of an unfortunate choice of coordinates, and system (2) is a 
canonical form for system (1). 

If we want to build a controller for system (1), then a controller can be 
designed for system (2) and applied to system (1) through the transformation (from 
x-space to y-space) and its inverse. This is the basic design technique for the 
automatic flight control of aircraft applied in references 2-5. The general multi- 
input version of transformations of nonlinear to linear systems is given in 
reference 6 . 

Suppose we retain the assumption that g,[f,g]» . . . , (ad n 1 f ,g) span IR n for 
points near the origin, but remove the assumption about involutivity . What canonical 
form can we then derive for system (1)? We certainly want this canonical form to 
show "a certain amount of linearity." Under what Lie bracket conditions does this 
canonical form exist, and can we actually present the form in such a way that it 
accentuates the important Lie brackets? 

For example, the system 


(3) 


can be transformed by y x = x x , y 2 = x 2 , y 3 = x 3 + x 2 , v,= u + 3x 2 (x 3 + x 2 ) to 

y 1 = y 2 + (y 3 - y 2 ) 2 = y 2 + y 3 - 2 y\y 3 +. y\ 

y 2 = y 3 

y 3 = v 

Letting y 1 = y x + (y*)/2, y 2 = y 2 , y 3 = y 3 , v = v, we have (deleting the ~ notation) 

yi = y 2 + y 3 + y 2 ' 
y 2 = y 3 
y 3 - v 

If we write system (5) as y = f(y) + vg(y) , system (5) cannot be reduced to the 
linear system (2) since the set {g, [f,g]> is not involutive (see ref. 1). However, 
because 




*1 

= k 2 

+ 

2 

X 3 

X 2 

” X 3 

+ 

x 2 , f = 

# 




x 3 

= U 





2 



( 6 ) 



~0" 


" 2 y 3 ~ 


~- 2 ~ 


"I + 6y 

g - 

0 

, [f,g] » 

i 

» [g» Ef >§]] = 

0 

, (ad 2 f,g) = 

0 


_1_ 




_0 _ 


0 


we can write the form (5) as 


~y~ 


"o' 


“<f 


“yf 


~y~ 


"Vf 

h 

= 

0 

+ 

y 3 

+ 

0 


0 

+ 

0 

i — 

v;. 

C*J 

L 


_v_ 


_o_ 


_0 _ 


_0- 


_o _ 


= gv + [f,g] 


=o 



(7) 


+ (ad 2 f , g) 


y 2 =y 3 =° ^2 


+ i(-(ad 5 [f,g],(ad 2 f,g)) 


y 2 =y 3 =o 72 


In equation (5) we point out the two dimensional linear subsystem which is a string 
of integrators 


y 2 = y 3 

y 3 = v 


( 8 ) 


Equation (5) (or (7)) is the canonical form of interest to use for this example. 
We derive such forms for n dimensional systems, and give conditions under which a 
form shows a particular type of linear subsystem (as in eq. (8)). We provide a con- 
structive procedure to move a system from the original coordinates to the canonical 
form, and this procedure is the main contribution of this paper. Also, the canonical 
form is unique only up to solutions of partial differential equations we introduce. 


These canonical forms are important because they show the intrinsic nonlineari- 
ties of a system and the obstructions to having a transformation to the linear sys- 
tem (2). Also, as in all areas of mathematics, it is often much easier to prove 
theorems if we assume some canonical form which exhibits the basic mathematical prop- 
erties (e.g., Jordan form, rational canonical form, companion form, Brunovsky (ref. 7) 
form) . 

As noted before, one interesting problem is characterizing those nonlinear sys- 
tems which can be transformed to controllable linear systems in canonical form (this 
research contains multi-input as well as single-input results). We refer here to the 
work of Krener (ref. 8), Brockett (ref. 9), Jakubcyzk and Respondek (ref. 10), 

Hermann (ref. 11), and the authors (refs. 1,6, and 12-15). Hermann applies the 
theory of equivalence of exterior differential systems to study canonical forms under 
feedback for nonlinear systems. 


3 


The results of this paper can easily be generalized to multi-input systems. We 
consider this problem and examine a two- input control system as an example at the 
end of the paper. 


II. DEFINITIONS AND PRELIMINARIES 

For vector fields f and g on IR n we define the Lie bracket of f and g 

[f.sl -ff 8 -ff f ' 

where 9f/9x and 9g/9x denote Jacobian matrices. This is the negative of the stan- 
dard definition, but it allows for easier notation in our present study. We can also 
define [f,[f,g]], [g,[f,g]], [f , [f , [f ,g] ] ] , [g, [f , [f ,g] ] ] , etc. In fact, we let 

(ad°f,g) = g 

(ad x f ,g) = [f,g] 

(ad 2 f,g) = [f,[f,g]] 


(ad k f,g) = [f , (ad k-1 f , g) ] 

A set of #°° (smooth) vector fields {X 1 ,X 2 X r } on JR n is involutive if 

there exist © functions so that 

r 

[Xi.X^Cx) = £ Yii k (x)X k (x) , 1 < i, j < r, i i j 

k=i 

If f is a vector field on IR n and h(x) is a ^ function, then 


L f (h) - <dh.f> t. . • 


where dh is the gradient of h. Similarly, for a #°° one form 
a) = Wj dx x + to 2 dx 2 + . . . + to^ dx n we define 

(co,f) = u^fj + oo 2 f 2 + . . . + to n f n 


and the Lie derivative of to with respect to f 




where * denotes the transpose and 9co*/9x and 9f/9x are Jacobian matrices. For g 
a vector field on IR n and f, h, and to as before we have the formulas (see ref. 12) 
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<«*), [f ,g]> 


Lf(w,g> = <L f (w),g> - 


(9) 


dLj (h) «* Lf (dh) 


( 10 ) 


The following lemma is implied by results in reference 1, but for the sake of 
completeness, we present a proof. Let f and g be as in system (1) (the assumption 
of real analytic can be relaxed) . 


Lemma 1 - Suppose the set of vector fields .{g, [f,g] , . . . , (ad n_1 f ,g)} are linearly 
independent and {g,[f,g], . . . ,(ad9f,g)} is involutive for some integer q, 

1 ^ q ^ n - 2. Then the sets of vector fields 


• •; . , (ad9 -1 f jg )} 

(g»[f»g]» ♦ • . , (ad9~ 2 f , g )} 


(g»[f.g]} 

are also involutive. 

Proof. Since (g,[f,g], . . .,(ad9f,g)} is involutive, by the classical Frobenius 
Theorem there are functions T X ,T 2 , . . • .T^q-i so that dT-^dT;,, . . ..dTn.,,.! 
are linearly independent and <dT i , (ad)f ,g)> = 0, i = 1,2, . . . ,n - q - 1 and 
j = 0,1, . . . ,q. There is at least one T^, say T n _ q _i, so that 

<dT n _ q _ 1 ,(ad9 +1 f,g)> ^0 


Define T n _ q - Lf(T n _ q _ 1 ) 


By formula (9) 


L f <dT n-q-i’( adj “ lf ’g>> = (dT n _ q _! ) , (ad^ _1 f ,g)> - (dT^^ , (ad^f ,g)> 

for 3 - 1,2, . . . ,q 
Thus 

< L f( dT n-q-l),( adj_1 f,g)> - 0 

and by equation (10) 

<dT n-q» ( ad ^" lf »g)> = 0 , j = 1,2, . . .,q 

Hence 

<dTi, (ad^f ,g)> = 0 , i - 1,2, . . . ,n - q and j = 0,1, . . . ,q - 1 


5 



Formulas (9) and (10) are applied again to show 

<dT n _ q ,(ad<lf,g» = <dT n _ q _ 1 ,(ad < l +1 f,g)> * 0 

We claim that the vectors dT^dT^, . . • > dT n -q are linearly independent, know- 
ing that dTjjdTg, . . •»dT n _ q _ 1 are independent by assumption. Take constants 
Cj , c 2 , ♦ • « ) c n _q so that 

Ci dT x + c 2 dT 2 + • • , + c n _ q dT n _q =0 
We dual product this with (ad^f ,g) 

c 1 (dT 1 , (ad^f ,g)> + c 2 (dT 2 , (ad9f ,g)> + . . . + c n _ q <dT n _ q , (ad^f ,g)> =0 

and find that c n _q = 0. Hence Cj = c 2 «■ . . . = c^-q-i = 0 and the desired gra- 
dients are linearly independent. 

By the Frobenius Theorem, linearly independent gradients satisfying 
(dT^, (adJf ,g)> = 0 , i = 1,2, . . . ,n - q and j = 0,1, . . . ,q - 1 

imply that the set {g,[f,g], . . . , (ad9 -1 f ,g) } is involutive. Repeating this process 
q - 2 more times completes the proof. 

A #°° distribution A on IR n is an assignment A (x) of a linear subspace of 
H n at each point x of ]R n . We assume that A is of positive constant dimension 
k and identify A with the set of vector fields in it. We also let A be involu- 
tive and regular (in ref. 16 this means that the quotient set IR n /A is a #°° mani- 
fold). The distribution A is (f,g) invariant if there exist a(x) and B(x) so that 

[f,A] C A 

[g,A] C A 

where f = f + got and g = g£. Here f and g are as in equation (1). 

The following lemma is found in reference 16. 

Lemma 2 . Let A be an (f,g) invariant regular distribution. Then A induces a 
regular equivalence relation on ]R n such that the dynamics in system (1) passes to 
the quotient denoted by lR n /A, whose dimension is n - k. 

In our later application of this result, since our theory is local, IR n is 
replaced by an open neighborhood of the origin. 


III. CANONICAL FORMS 


We consider system (1) and assume the following conditions hold on an open 
neighborhood of the origin in ]R n . The vector fields g,[f,g], . . . , (ad n_:L f ,g) are 
linearly independent and the set {g,[f,g], . . ..(ad^g)} is involutive for some 
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fixed integer q, 0 < q < n - 2. All arguments and results are local and hold in a 
neighborhood of the origin. Some coordinate changes used are similar to those in 
reference 9. 

Since g is nonvanishing there is a well known coordinate change so that g 
becomes 

r:i 


We thus assume that equation (1) is 



Using feedback (new u = old u + f n ) we can assume f n = 0, 


(ID 


Let 


x 


■ktbu.JP 


' 3x 


X + gu 


( 12 ) 


be the linearization of system (11) about the origin. Since at the origin 
g = b,[f,g] = Ab, . . . ,(ad n-1 f,g) = A n-1 b, we have that this linearization is con- 
trollable. It is well known that coordinate changes and feedback can be applied to 
take system (12) to the string of integrators (2). Hence equation (11) can be put 
in the form 


rvi 


m m 

x 2 + . . . 


mm ■» 

0 

CM 


x 3 + . . . 


0 

• 


• 


• 

• 

= 

• 

+ u 

• 

• 


• 


• 

*n-i 


Xjj f . . . 


0 

A. 


0 


-1- 


f + ug 


(13) 


where + . . . denotes higher order terms. 
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Replace ^ + . . . by a new x^ compute Xjj, and apply feedback to return, the 
last entry in f to 0. For equation (13) we have 




x 2 + . . . 


V 

*2 


■ Xg * 


0 

• 


• 


• 

• 


• 


• 

• 

= 

• 

+ u 

• 

*h-2 


x n-i + • • • 


0 

• 





*n-l 


% 


0 



0 


1 


= f + ug 


(14) 


We consider the terms 


a 1 (.x 1 ,x 2 , . . . ,x n _ 1 )x n 
^2 ^1 * * • • »X n -l)x n 

a n-2^ x i» x 2> • • *> x n-i) x n 
0 


in the f vector field. Using coordinate changes on x 1 ,x 2 , 
take 


a i( x i> x 2> • • • #x n-i^ x n 


"o 

a 2 ( x i 9 X 2 9 • • ,x n*-i) x n 


0 

• 


• 

• 

to 

• 

• 

a n _ 2 (x 1 ,x 2 , . . • > x n _ 1 )x n 


0 

x n 


1 

_ 0 


0 




space 


we 


Thus we have form (14) where a x = a 2 = . . . = a n _ 2 = 0. 

We have not used the fact that the set {g,[f,g], . . . ^(ad^g)} is involutive 
if q > 0 and the implications of Lemma 1. We take (14) and compute certain Lie 
brackets. 
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We find 


0 


9x n 


I 

<M 

<*> 

1 

0 


9f 2 


9 2 f 2 


9Xn 


9x 2 

• 

• 

. [f.gl = 

• 

• 

» [[f »g] *g] - 

• 

• 

' • 


• 


• 

o 


9f n -i 


3 2 fn-i 



9x n 


9x n 

l 


0 


° 


with 9f n _ 1 /9x n = 1 and 8 2 f n _ 1 /9xn = 0. By Lemma 1 we must have Yx (x) and y 2 (x) 
such that 

[tf.gl.g] =Y 2 (x)[f,g] + Yi(x)g 

This implies that [[f»g],g] = 0 vector and 

9 2 f 1 /9x£ = 9 2 f 2 /9xn = • . . = 9 2 f n _ 1 /9xn = 0 (and hence 

9f x /9x n = 0, 9f 2 /9x n = 0, . . 9f n _ 2 /9x n = 0 since a 1 = a 2 = . . . = a n _ 2 = 0). 


Replace x n-1 + . . .in (14) by a new x n _ 1 , compute > make coordinate 

changes for x n , compute x n , and apply feedback to return the last entry in f to 0. 
For equation (14) we have 


x i 


^2 « 


'o' 

x 2 


x 3 + . . • 


0 

• 


• 


• 

• 


• . 


• 

• 


• 


• 

CO 

1 

•*F 

— 

x n-2 • • • 

+ u 

0 

x n -2 


x n _i 


0 

Xn-i 


x n 


0 

1 

i 


0 


1 

™ m 


= f + ug 


(15) 
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We examine the terms 


b x (xi ,x 2 , . . • > x n~2^ x n-’i 
b 2 (x 1 ,x 2 , . • • »x n _ 2 )x n - 1 

• . % 

^n-3^ x i» x 2* * * ** x n- 2 ^ x n-i 

*n-i. 

0 

0 

in the f vector field from equation (15). Coordinate changes on 
Xi ,x 2 ) « • • jX ^|_2 space move 


b 1 (xi,x 2 , . . . ,x n _ 2 ) 


"o' 

b 2 (x^,x 2 , * . • 


0 

• 


• 

• 


• 

• 


• 


to 


^n-3^ x i» x 2* * * • » x n-2^ 


0 

1 


1 

0 


0 

0 


_ 0 


Hence in equation (15) we assume b, = b, = . . . = b , = 0. 

x £. n- 3 

^ *1 > then { g» [ f , g] » (ad 2 f , g) } is involutive by Lemma 1. Computing again. 


0 

; 

. 

0 


> £ i 


"Vf-f 




3x n _i 



0 


0 




3 2 f 2 



^n-i 


K-i 

4 

• 

» [f.g] = 

• 

• 

, (ad 2 f,g) = 

• 

• 

[(ad 2 f,g),(f,g]] = 

• 

• 

0 


• 

0 


• 

1 


♦ 

0 

0 


1 


0 


0 

1 


0 

| 

0 


0 
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There are y l (x), y 2 (x) , and y 3 (x) such that 

[(ad 2 f,g),[f,g]] - Y 8 (x)(ad 2 £,g).+ Y 2 (x)[f,g] +Y 1 (x)g 

But this implies that 3 2 fi/9x^_ 1 = 0 for i - 1,2, . . . ,n - 3, and hence 

3f i /3x n-i = °» 3f 2 / 9x n-i " °» . • • •» 9f n -3 /9x n-i = 0 since 
b x » b 2 = . . . b n _ 3 = 0. 

We then repeat the above arguments until our equation (15) becomes 


*1 


" f i “ 


mm wm 

0 

X 2 


fa 


0 

• 


• 


• 

• 


• 


• 

• 


• 


• 

CM 

1 

O* 

1 


f n-q -2 

+ u 

0 

^-q-l 


^n-q-i 


0 

• 


• 


• 

• 


• 


• 

» 


• 


• 

> 

1 

M 


f n-x 


0 

1 

1_ 


0 


1 


with 


00 

X 2 + E 

i=2 

h i,i.i x i + 

i U] 11 

00 

E 

i=i 

h i,j,i^ x i» x 2» * 

. . ,Xj_ 1 )xj- 

. 00 

x 3 + E 

h 2 , 1 , i x i + 

n-q 

E 

00 

E 

^2 , j ,i ^ x i ,x 2 5 • 

• •» x j_ 1 )xj 

±=2 


j=2 

i-i 



“ i n-q oo 

^n-q -2 “ x n-q-i + E h(n-q- 2 ) , 1 , i x i + E E ^(n-q- 2 ) , j ,i ^ x i » x 2 > • • •> x -j_i) x ^ 

1=2 j=2 1=1 J J 

f n-q-i = ^n-q 


^n-i x n 


f n " 0 


(17) 


the definitions and functional dependences of the h's being obvious. 
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We apply coordinate changes on x lf x 2 , . V . ,x n _ q _ 1 space to take 


h i,n-q,i^ x i» x 2» • • • » x n _q_i) 


“cf 

^2,n-q,i^ x i»- x 2» • • f» x n -q-i) 

■ 

0 

• 


• 

• 


• 

• 


« 

h (n-q- 2 ),n-q,i( x i> x 2 » • • • » x n -q-i) 

to 

0 

1 


1 

0 


0 

• 


• 

• 

• 


• 

1 

o 

1 


J> m 


Also in new coordinates on x ,x , . 

h i,n-q-i,i< x i’ x 2’ * 

h 2,n-q-i,i^ x i‘ x 2» * 

' • 

• 

h(n-q-3) ,n-q-i , 1 ( X x » x 2 » 
1 + h (n-q-2) ,n-q-i,i( x i» x 

0 . 

0 


• » x n -q -2 space 


• 

** 

X 

t 

1 

PO 

'w' * 

1 


”o~ 

* > x n-q~2^ 


0 



• 

« 

• ■' 

• • • > x n-q-2^ 

appears as 

0 

* * * “ x n-q-2 ) 


1 



0 

• 



• 

• 



0 


Similar arguments allow us to arrive at 
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where 



( 18 ) 


1 = X 2 + 2 X) ^1 , j , i ( X 1 » x 2 » * * •> x i_i) x '4 
j=l 1=2 


f 2 = X 3 


w n*"Cj 00 

+ S ^2,l,i x l + S ^2 2 ,1 ^ x l^ x 2 + iC ^2 , j ,i ( X 1 » x 2 > * • •> x j-l) x l 

1=2 1=1 ' * 1=3 1=2 J 


n-q-2 oo 


n-q-2 


^n-q-l + 11 (n-q-2 ) ,i,i x i + h (n-q-2) , j ,1 ( x i > x 2» • • • > x j-i) x j 

n-q oo 

+ 53 23 h (n-q- 2 ) , j ,i^ x i ,x 2 ’ * * ,,x 1-i^ x i 

j=n-q-i i =2 j 


n-q-i 


f n-i " x n 


f n = 0 


with new h functions. 

Equation (18) is our canonical form for the nonlinear system (1), 
the following result. 


We have proved 


13 



Theorem 1 . Suppose f and g in system (1) satisfy, in a neighborhood of the origin, 

g»[f»g]» • • . ,(ad n-1 f,g) are linearly independent, 

2- g.[f,g] (ad^f ,g) are involutive. 

Then there are new state space coordinates, which are also called x^Xj,, . . . ,x_, 
and a new control u so that system (1) becomes system (18). 

If q = n - 2 then equation (18) is a linear system, and we have the results of 
reference 1. 

Suppose we are given two systems, (18) and 


~ *1 ~ 


- ^ 


“o’ 

y 2 


^2 


0 

• 


• 


• 

• 


• 


• 

• 


• 


• 

• 

yn-q-2 



f 

^n-q-2 

+ V 

0 

^n-q-i 


^n-q-i 


0 

• 


• 


♦ 

♦ 


• 


• 

y n -i 


^n-i 


0 

__ y n _ 


0 


0 


with f 1 »f 2> . . • »f n _! as in equation (19)_except that x 1 ,x 2 , . , ..x^u are 
replaced by yi>y 2 > • • • >yn> v an d b by h. We wish to know if they are equivalent 
in the sense that system (20) is simply system (18) in new coordinates 
yi (x) ,y 2 (x) , . . . ,y n (x),v(x,u). If this is true, there is a transformation from 
system (18) to system (20). Since involutivity is invariant under these transforma- 
tions we assume that the integer q is the same for both systems. 

If (18) is equivalent to (20) , then we have a transformation as above which 
satisfies the following partial differential equations. 
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^ 9yx 

£ £i ■ £ ‘ 


A 9y 2 

£ £± ' £i 


n dy , 

V' ^n-l . 7 

axi fi " f n-i 


( 21 ) 


Here we need y lf y 2 y n -i as functions of x 1 ,x 2 , . . . .x^ only and y a 

function of x 1 ,x 2 , . . .,3^. Writing these equations out we obtain 

3Yi 3y x 3y x 

3x * * • ) + 0 X ( x 3 + • • • + Xt 


3x ~ n 
dX n-i 


n-q 00 

= y 2 + E £ h l,j,i (y i‘ y 2 y i-l )y j 

j=l i=2 J J 


3y 2 

9xi 


(x 9 + 


= y : 


9y 2 3y 

* ^ ^ 9x7 ^ X 3 * * • ) + • • • + TTT x. 


3x n 

ox n-i 

n-q oo 

j- 


” _ ” . n-q oo 

+ £ h 2,i,i y i + £ fi 2 ,2,i(yi)y2 + E E S 2,j,i(yi.y 2 . • • 

i=2 1=1 . -1 = 3 1=2 j J 


9y l 3y 

(x, + . . .) + ~~ (x 3 + . . .) + 


9x n 


9x, 


9y n-i 

9xn_i 


*n 


r n 


( 22 ) 


In the left-hand side of the first (n - 2) equations, the terras 
(3y 1 /9x n _ 1 )x n , (9y 2 /9x n _ 1 )x n , . . . , (3y n _ 2 /9x n _i)x n consist of functions of 
x i > x 2 » • • • » x n-i times the variable x,^. In the right-hand side of these equations, 
the only possible terms of this type are those that contain a y n (and this can occur 
only if q = 0). But such terms are raised to powers of two or greater. Hence 
yi»y 2 » • • • *yn -2 are functions of x 1 ,x 2 , . . . .x^j only. 

£ We examine the terms O yi /dx n _ 2 )x n . x , (9y 2 /9x n . 2 )x n _ 1 , . . . , (9y n > 3 /9x n _ 2 )x n - 1 
from the left-hand side of the first (n — 3) equations* The only possible way such 
terms can appear on the right-hand side of these equations is through y n-1 and y n 
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(if q 5 1), but again these are raised to powers greater than one. We repeat this 
argument to achieve 

y 1 (x 1 ,x 2 , . . .,x n _q_ 2 ) 

y 2 * X 2 * * * • »Xn_q— g) 


y n-q-2^ x i» x 2’ * • • >x n-q-2^ 
yn-q-i^ x i» x 2» • • ,,x n-q-x^ 


(23) 


( x i » x 2 * * * *» x n-x^ 

Yn ( x i» x 2» • • ••Xq) 

The partial differential equations (22) now become 
n-q-2 9 7i 

E sr (*i+i + • • •) 

i=i i 


n-q » 

= y 2 + E E iI i > j,i<y 1 *y 2 » 

j=l i=2 




n ~q -2 3 y, 

E 9^7 (x i+i 

i=>x i 


w _ >» n-q 00 

= y 3 + E H 2 } 1 >i yj- + E S 2, 2 ,i(yi)y2 + E E E 2 ,j,i(yi.y 2 * • • •»yj-i>yj. 


.n-q- 2 3 V 

oy n-q-2 f 

E 9 xj ^ X i+i + * * *) 


i=l 


” n-q-2 oo 

y n-q-x + E S (n-q-2 )fl ,i7i + E I) Nn-q-a) ,j .i^x * * * »yj- 1 )yj 

n-q oo 

+ E E h(y 1 ,y 2 , . . . ,y i _ 1 )yf 

j=n-q-x i=2 J 


(24) 
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£ 


3x^ 


vx i+1 


ii-q-l 


n-q-2 3y 

£ -si* <*i« + 


• • -) + . ^ , ~s, 

i=n-q-i i 


n ~q 3y« „ 

^ n~q _ 

ix^ x i+i ^n-q+i 


(24) 

(Cont) 


n -f 2 9 y n -i 

i=i 


ox,- 


( x i+i + . . .) 


n-l 

+ E 

i=n-q-i 


^n-i 

3x^ 


*1+1 



Thus we have reduced the equivalence problem of two systems to that of finding a 
transformation of the form (23) which satisfies equations (24). We present examples 
in this direction. 


Example 1 . We take two systems which are already in the canonical form. 

* . 1 2 
X 1 = X 2 + x l x 3 

x 2 = x 3 
X* = u 

h = y 2 + y^r 

y 2 = y 3 
y 3 = v 


(25) 


(26) 


Now by the results of reference 8, there is no state space coordinate changes 
which takes one system to the other. However, if we allow transformations also 
involving controls, this is possible. 


By equation (23) we have y^x-^, y 2 (x 1 ,x 2 ) and y 3 (x, ,x 2 ,x 3 ) . Substituting into 
equation (24) we find 


Since we are looking for 

y 2 ( x i » x 2 ) = 

y 3 (x 1 ,x 2 ,x 3 ) = 



ya + yiy 3 

y 3 


/ 


real-analytic solutions we expand 

r 2o( x l) + r 2l( x i) x 2 + r 2 2 ( X 1 ) X 2 + * * * 

r 30 ( X l » X 2 )' + r 3 l < X i ’ X 2 ) X 3 + r 32 ( X X » X 2 ) X 3 + • 


(27) 


(28) 
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Thus the first equation in (27) gives 


fox _9y 

9x x x 2 + 8x 


9yi (i \ 

8x x \2 X l X 3y = r 20^ X l^ + r 21^ X l^ X 2 + r 22^ X l^ X 2 + • • • 

+ yi< X l>( k r 30< X 1 ’ X 2> + r 3 1 ( X 1 ’ X 2 > X 3 + r 32 (X l» X 2 )X 3 + -'•)* 


Computing we have 

y2< x i» x 2> " r 2 1 ( X X ) X 2 1 
y 3 ( x i» x 2 » x 3 ) " r 31 (x 1 )x g 

and from equation (27) 

9y i 1 2 / x : 

g x X 2 + g x 2 XlXs " r 21^ X l^ X 2 + yi< x l)( r 31< x l>*3) 


9y 2 7 ,1 2 \ . 9y 2 , . 

g^- ^X 2 + 2 x l x 3 j + 3^- x 3 - r 3x (x 1 )x 


3x x 2 

9y 2 


Hence 9y 2 /8x 1 - 0, y a (xj,x 2 ) = y 2 (x 2 ) - r 21 x 2 , 3y 2 /9x 2 ==-r 31 (x!) - r 31 , and 
3y 1 /9x 1 = r 21 where r 21 and r 31 are constants. Integrating we have 

y x = r 2X x x 

y 2 = r 31 x 2 

y 3 - r 3 x x 3 

Since y 2 = r 31 x 2 , r 3l = r 21 and there is a constant r / 0 such that 

y x " rx X 


y 2 


= rx. 


y 3 = 


= rx. 


To find r we substitute into the first equation in (30). 


or 


rx 2 + rx 1 x 3 = rx 2 + rx x r 2 x 2 


r = r' 


( 29 ) 


(30) 


with solution r = l/(2) a / 2 . 
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Hence the transformation 


Yl ~ Si * l 


y 3 = 


/2 


V = — U 

Si 

takes system (25) to system (26). 

Example 2 . Again our systems are in canonical form 


X 1 

= X 

• ■ 

X 2 

= X 

• 


X 3 

= u 

• 

yi 

= y 


C 2 X 3 


■- *2 T y x y 2 x 

y 2 = y 3 
y 3 - v 

The partial differential equations (24) in this example are 

9y x 


9x n 


(x 2 + x 2 x 3 ) = y 2 + Y x y\ 


9 y 2 , 2x 9y 2 

jsr <x * + vJ + 


* 3 = y 3 


(31) 


(32) 


(33) 


We take the expansions (28) for y 2 (x is x 2 ) and y 3 (x x ,x 2 ,x 3 ) . The first equation in 
(33) is 

3y 3y x 

9^7 X 2 + 3^ X 2 X 3 = r 20< X l> + r 21<V X 2 + r 22 (X ! )X 2 + * * * 

yi(Xi)(r 30 (x 1 .X 2 ) + r 3 i ( x i » x 2 ) x 3 r 3 2( X l» X 2^ X 3 '*'•••) 
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This implies 


y 2 (x 1 ,x 2 ) - r 21 (x x )x 2 
y 3 (Xi»x 2 ,x 3 ) = r 31 (x 1 ,x 2 )x 3 


( 34 ) 


and from equation (33) 

3yi 


3x, + 3^ X 2 X 3 = r 2i ( x 3,) x 2 + yi(x 1 )r| 1 (x 1 ,x 2 )x| 


3y 2 3y 2 

toT < X 2 + X 2 X 3> + 3lT X 3 " r 3l(x X >X 2 )x 3 


(35) 


Therefore, 3y 1 /3x J = r 21 (x 1 ), (3y 1 /3x 1 )x 2 = y^^x x^, 3 y 2 / 9x ! ■ °» y 2 ( x 2 ) 
3y 2 /3x 2 = r 31 (x 2 ), where r 21 is now a constant. Also 

r 2 i(xi)x 2 = y x r 3 i(x 2 )] 
r 2i - r 31 < x 2> 

Then r 31 is a constant and r 31 = r 2 ,. The first equation in (36) 

r 21 X 2 = r 3xyi< X l> 


r 21 X 2’ 


(36) 


is impossible to satisfy, and our systems are not equivalent. 

We now discuss the linear subsystem contained in equations (18) and (19) . It is 
clear that {g,[f,g], • . . ,(ad^f,g)} being involutive implies the existence of the 
linear subsystem 


x n-q-i “ ^-q 
x n-q = x n— q-i 

(37) 

Vi " *n 
in = u 


However, there may be a "larger" linear subsystem as illustrated by the following 
example. 
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Example 3 , Consider on ]R M 


x, = x„ 


x„ 


x lt = u 


+*i\ 


(38) 


Now the set {g,[f,g]} is not involutive, implying q = 0. But the coordinate changes 




y 2 = x 2 

y 3 - x 3 +x a 

y 4 = + 2x 2 x 3 + 2x 2 

v = u + 2x 2 x 4 + 2(x 3 + x 2 )x 3 + 6x 2 (x 3 + x 2 ) 


yield the system 


*1 = y 2 + yl 
y 2 - y 3 

y 3 = >4 

y 4 = v 

Thus we have a ‘'larger 11 linear system than indicated by the involutive assump- 
tion. However, the coordinate changes just used are contained in those applied in 
the proof of Theorem 1. Hence the process we introduced in Theorem 1 will provide 
the linearity beyond that given to us by the integer q. In our system (39) above we 
can also use coordinate changes on y 1 »y 2 >y 3 space to send 


l4 y2y 3 " 


"o' 

0 

to 

0 

i 


1 

0 


0 

m « 


4y 2 y 3 y 4 + 4y§yJ 


(39) 


For system (39) let A be the distribution consisting of the vector field 


1 

0 

0 

0 
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Then [f,A] C A and [g»A] c A and A is an (f,g) invariant distribution with a = 0 
and 8 = 1. By Lemma 2, the dynamics on the quotient manifold IrVa is linear. 

We return to system (18) and let A be the distribution spanned by the vector 
fields 


“l“ 


“o" 


“ 0 

0 


1 


• 

• 

0 


0 


• 

0 


9 


9 • • • 

1 

• 


• 


0 

• 


• 


• 

• 


• 


• 



_ 0_ 


• 

_ 0_ 


where the 1 in the last vector field is in the (n - q - 2) position. Since the 
Jacobian matrix of f in equation (18) has only zero elements on and below the 
diagonal, 

[f.A] C A 

Trivially, 

[g»A] c A 

and A is an (f,g) invariant distribution with a = 0 and 8=1. 

Hence the dynamics on the quotient manifold ]R n /A, where lR n is actually an 
open neighborhood of the origin, is linear. This dynamics is actually the linear 
of equations (18) and (19) , and linear design techniques can be applied on 
TR /A. As Example 3 shows, if more linear dynamics exists, this procedure can still 
be applied. 

Equations (18) and (19) emphasize the essential nonlinearities of the nonlinear 
system (1). Since a part of the system actually appears in linear form, a design 
technique based on equation (18) should be simpler than on the original system. 

If the integer q is equal to n — 2, then equation (18) is: a linear system, 
and we have a transformation from a nonlinear to a linear system. The design of an 
automatic flight controller involving a multi-input generalization of this transfor- 
mation is presented in references 2-5, as mentioned before. 

Suppose we have a multi-input system 

m 

£(t) = f(x(t)) + £ u i (t)g 1 (x(t)) (40) 


where fjg^ . . •jgm are real analytic. It is not difficult to extend the results 
of this paper to develop a canonical form for equation (40). We consider the case 
n = 7, m = 2, and g x , [f >g x ] , (ad 2 f ,'gj) , (ad 3 f ,g x ) ,g 2 , [f ,g 2 ] , (ad 2 f ,g 2 ) span ]R 7 . 
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We also assume the set {gi ,g 2 , [f ,g x ] , [f ,g 2 ] > is involutive, and a parallel result to 
Lemma 1 implies the set {gi»g 2 } is involutive. Thus we take (possibly renaming 
controls) 


Si = 


0 


"o 

0 


0 

0 


0 

1 

and g 2 = 

0 

0 


0 

0 


0 

0 


1 


After applying feedback we have as our system 


• 

x i 


U-i 

1 


“ “1 
0 


"o’ 

x 2 


fa 


0 


0 

*3 


u 


0 


0 

K 

« 

0 

+ u x 

1 

+ u 2 

0 



f 5 


0 


0 

*6 


f 6 


0 


0 

js_ 


0 


o 


1 


(41) 


We use linear feedback and linear coordinate changes on ]R 7 to take the lineariza- 
tion of (41) about the origin to Brunovsky (ref. 7) form. Hence equation (41) 
becomes 


*1 


x 2 + . 

• • 


V 

X 2 


x 3 + . 

• » 


0 

*3 


x 4 + . 

• • 


0 


- 

0 


+ u, 

1 

X 5 


x 6 + * 

• • 


0 

*6 


X ? + . 

• * 


0 

3_ 


0 



0 


+ u, 


(42) 


where + . . . denotes higher order terms. 

Replace x 7 + . . . by x 7 , x 4 + . . , 
feedback so that equation (42) becomes 


by x 4 » compute x 7 and x 4 , and apply 
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( 43 ) 
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The set {g x ,g 2 , [f,g 1 ],[f,g 2 ]} being involutive implies that S^/Sx 2 , 3 2 fi/3x^, 
3 2 fi/3x 4 3x 7 vanish for i = 1,2,5. Also, coordinate changes can be made on 
Xj^Xjj.Xg.Xg.Xg space to convert 


[f,g x 3 to 


V 


*0" 

0 


0 

i 


0 

0 

and [f,g 2 ] to 

0 

0 


0 

0 


1 

.0. 


c m 


We assume (43) is in these new coordinates. 


by a new x 6 , compute x 3 and x 6 , 


Replace x 3 + . . . by a new x 3 , x 6 + . . 
make coordinate changes for x 4 and x ? , compute x^ and x ? , and apply feedback to 
return the fourth and last entries in f to zero. Thus we find 


with 


! 

i-4 

• X 

1 


V 


0 


V 

X 2 


^2 


0 


0 

X 3 


f 3 


0 


0 

x 4 

= 

U 

+ 

1 

+ u 2 

0 

X 5 


f 5 


0 


0 

X 6 




0 

! 

0 

x 7 


^7 


0 


1 


f i “ x 2 + ^1 


f 2 = x 3 


t 3 = x 4 
f 4 = 0 
fs " X 6 


f 6 = X 7 


f 7 = 0 


(44) 


(45) 


f x being a function of x 1 ,x 2 ,x 3 ,x 5 ,x 6 only and containing no linear terms. 
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It is interesting to note the various equations (e.g., eq. (44)) that result if 
we take different involutive assumptions (in the above we took {gi,g 2 , [f ,gi] , [f ,g 2 ] } 
as our involutive set). For example, if the sets {g ; ,g 2 . [f .gj.] , [f ,g 2 ] ) and 

{gi»g 2 > [f »gil » [f »g 2 ] » (ad 2 f ,g x ), (ad 2 f ,g 2 )} are both involutive, we have a linear sys- 
tem as shown in reference 6. 


IV. CONCLUDING REMARKS 


We have introduced a canonical form for the nonlinear system 

x(t) = f(x(t)) + u(t)g(x(t)) 

which emphasizes the intrinsic nonlinearities of the system. This form is derived by 
proceeding through a series of coordinate changes in state and control Variables. 
Applications of this canonical form theory to the problem of system equivalence and 
an extension to multi-input systems are also discussed. If a system is transformable 
to a controllable linear system, then its canonical form is the Brunovsky one. 
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